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We describe a method to perform a single quantum mea- 
surement of an arbitrary observable of a single ion moving in 
a harmonic potential. We illustrate the measurement proce- 
dure with explicit examples, namely the position and phase 
observables. 

PACS: 03.65.Bz, 42.50.Wm 



According to Quantum Mechanics every observable 
A is represented by a Hermitian operator A, and the 
eigenvalues ak of this operator represent possible out- 
comes of a single measurement of ^ Q . Given a state 
10) of our system the probability for obtaining is 
-Pfc = KV'fcl'/')^ where jV^fc) is the eigenvector associated 
with flfc, AjV'fc) = CLkli^k) While Quantum Mechanics 
states that any observable can be measured in princi- 
ple, in practice for a given system only a few observables 
like momenta, energy etc. are accessible in the labora- 
tory. In the present letter we will show how in a model 
system of an ion moving in a harmonic trap the mea- 
surement of an arbitrary observable of the ion motion 
can be implemented. In this way one can perform single 
measurements of observables thus far considered inacces- 
sible, such as the Pegg-Barnett phase operator , or any 
combination of position and momentum, such as angu- 
lar momentum etc. We emphasize that in contrast to 
quantum tomography 0|, where the full density matrix 
of a system is determined by repeated preparations and 
measurements, we are considering a single measurement 
of an arbitrary observable. 

Trapped ions are exceptionally well suited systems to 
study fundamental aspects of quantum mechanics ||^ . By 
coupling laser light to the internal degrees of freedom, a 
single ion can be cooled to the vibrational ground state of 
the trap, and its motion can be manipulated coherently 
to generate nonclassical states of motion. State mea- 
surement of the internal degrees of freedom can be car- 
ried out with essentially 100% efficiency using quantum 
jump techniques Q. These unique properties have stim- 
ulated a series of fundamental experiments from genera- 
tion of nonclassical states of motion |^ and Schrodinger 
cat states [Q, to implementation of quantum gates [gj. 
In addition, there have been proposals for tomographic 
measurements of the atomic motional density matrix , 
and quantum reservoir engineering [l0|] . 

To carry out the single measurement of an arbitrary 
observable A we will employ the following tools which. 



as will be shown below, are readily implemented in an ion 
trap: first, the ability to synthesize any motional state 
lip) of the ion starting from the vibrational ground state 
|0), ie. IV') = CA/;|0) where is a unitary operator; 
second, the ability to perform a filtering measurement 
to distinguish whether the ion is in the vibrational 
ground state or not. 

To illustrate the procedure we assume the ion to be 
in the unknown state pTp^ ]. Let us denote hy Uk 
the unitary time evolution operator which generates the 
eigenstate |-0fc) = Uk\0) of the observable that we want 
to measure. We first transform the state of the ion by 
the inverse transformation to |xo) = UqIc/)), and measure 
whether the ion is in the ground motional state or not, 
after which Uq is applied. If the ion was found to be in 
the ground motional state, the state after this procedure 
will be {ipo), which implies that we have measured ao- 
The probability for this to happen is 



mxo)\' = \{M<f>)\'^Po- 



(1) 



which agrees with the probability Pq of measuring . In 
the case that we do not measure the ion in |0), its state 
is projected onto 



|0o> = (l-|^o)(V'o|)|</')/|| 



(2) 



where || . . . || serves to normalize the state. We then trans- 
form the state of the ion according to |xi) = Ul\4>o), 
measure whether it is in the ground motional state or 
not, and apply the unitary evolution operator Ui. If the 
ion is found in |0), the state after this step will be 
The corresponding probability will be equal to 



(3) 



i.e. it coincides with the probability of measuring oi. In 
the case that we do not measure the ion in |0), we apply 
the unitary evolution operator C/j, and continue in the 
same vein. After k steps, the probability of measuring 
the ion in its ground state will be Pk, and the state of 
the ion will be projected onto the state ji/'fe)- In this way, 
we will measure the ion in one of the possible eigenstates 
\4^k), which is equivalent to obtaining as a result the value 
ak of the observable A. In practise this procedure must 
be carried out in a finite number of steps, so we restrict 
ourselves to a finite dimensional Hilbert space. 

Let us now show how the above procedure can be im- 
plemented in the case of an ion trapped in a harmonic 
potential. For the sake of simplicity, we will consider the 
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one-dimensional case; the generalization to more dimen- 
sions will be described elsewhere. Here we will discuss 
how to prepare arbitrary states of the ion motion out 
of the ground state |0), i.e. how to implement and de- 
rive the unitary evolution operators U using laser pulses. 
This method is conceptually similar to the one proposed 
by Law and Eberly in the context of cavity QED [^ . 
We will also simulate numerically the measurement pro- 
cedure described above for the case of a trapped ion. 

We consider a single two-level ion of ground level \g) 
and excited level |e), where spontaneous emission is negli- 
gible 1^ , trapped in a harmonic potential and interacting 
with laser fields. For a travelling wave configuration the 
Hamiltonian of this system in the interaction picture is: 



Htr — 



(4) 



where kx{t) = -q (de~"^* -I- dTe"^*), A = wj;, — wo is the 
detuning, is the laser frequency, ujq is the frequency 
gap between the ground and excited atomic levels, v is 
the trap frequency, 77 is the Lamb-Dicke parameter, fi is 
the laser Rabi frequency, and d^ and a are respectively 
the phonon creation and annihilation operators. 

The method for synthesizing states consists of starting 
out with a cooled ion in state jg, 0), and by an appropri- 
ate sequence of laser pulses to coherently distribute this 
amplitude to formed a desired arbitrary superposition: 
U\g,0) — Y^n=o^n,\9,n) ll^, where we are restricted to 
TLm+1- The calculational method assumes we are car- 
rying out the opposite; we begin with a given arbitrary 
superposition and coherently coalesce it into the state 
\g, 0). This is carried out as follows: given that the high- 
est value of n is N , we use a laser pulse to push all of the 
population of \g,n) into |e,rt — 1) (a "diagonal" pulse), 
and then all of the population of |e, rt — 1) into \g, n — 1) 
(a "vertical" pulse), and so on, keeping careful track of 
what is happening to all of the other populated levels of 
the ion all the while. Inverting the unitary operations de- 
scribing these pulses and applying them in reverse order 
will give us tj , as desired. 

For the "vertical" pulse, we tune the laser on reso- 
nance (A — 0). When is sufficiently small we can as- 
sume that there are no off-resonant transitions ||], that 
is to say the level structure of the ion can be consid- 
ered to be a series of parallel, isolated, two-level systems. 
Precisely what is meant by "sufhciently small" will be 
detailed later. Assuming this, and transforming to a ro- 
tating frame (t/rfv = e*""*"') one can easily show that the 
Hamiltonian (Q) reduces to: 
00 

^tr~Y.2 (""^^^"'^ + W){nl (5) 

n=0 

where rj„ = f7(n|e"'''("+°*) |n) is the effective Rabi fre- 
quency, given by: 



0„ = f^^ 



fc=0 



k\ ■ 



(6) 



The general unitary operation describing transfer of 
population between a ground and an excited atomic state 
(essentially a rotation) can be described by two parame- 
ters, 9 and Xj defined by: 



U = 



sm ( 



-e'-^ sin0 
cos 6 



(7) 



where \g) = (0,1), |e) = (1,0). For U{rge''^^\g) + 

TeC^'^'^le)) = y^r^ + r'^e^^f \g) (the desired operation in 

this case), tan 6* = Te/vg, and x = V'e — V'g- V'/ is the 
phase of the amplitude of the finally populated state. 

From (0) and (|) we can deduce that ^nt/'^ = 9W n > 
0. Thus, if we calculate 9 to clear the ground state of 
the "top" two level system (that system still populated 
having the largest number TV of phonons) 9^, for those 
systems possessing fewer (n) phonons, 0„ = ^,^9^ /^iq . 
For each two-level system (f) — x + 7i'/2. 

In the case of the "diagonal" pulse, we set A = —v 
(tune the laser to the lower sideband) and transform to 

a different rotating frame {Urid = e"'(°*°+''^/^)*). In this 
case, the Hamiltonian (^) can be approximated by 



00 „, 



1| +H.C. 



(8) 



n=0 



where VL'^ = VL{n\e-''^i^+^'')\n + 1) is the effective Rabi 
frequency, given by: 



.ai:(n;) ''"''^"''" <»' 



Note that the effective Rabi frequencies i7„ and fi^ are 
calculated from the full expansion of the exponential 
^-ikxt ^ and that our proposed scheme is therefore not re- 
stricted to the situation wheren <C 1 (Lamb-Dicke limit) , 
in contrast to previous work |^ and cavity-QED ||l3| . 
Referring again to (0), the desired operation is 

U{rge^'l'^\g) -t- ree"^=|e)) = J rj + rle^'l' ' \e) , so that 



tan0 = Tg/re, x — "^e ~ V'g + We thus deduce: 
= V 71 > 0; 6'n (apphes to the system of \g, n+1) 
and \e,n)) = il'^9N /fl'j^; 4> = X- To invert ^ for both 
pulses we simply add tt to the calculated x- 

Let us consider when the use of the approximate 
Hamiltonians which express the ion as a simple sum of 
non-interacting two-level systems is justified. For the 
sake of simplicity, we will give an estimate assuming the 
Lamb-Dicke limit, defined as 77 <C 1. In this case the con- 
dition ttotai^E ^ 1 must be fulfilled, where Ai? is the 
level shift induced by the next most important (off reso- 
nant) pulse when the laser is applied. Using second-order 
perturbation theory, this can be calculated to be F^/A, 
where V is the effective Rabi frequency of this secondary 
pulse and A is the detuning from the frequency gap of the 
two levels it couples together. The time taken for each 
"vertical" pulse t — 29/VL, and for each "diagonal" pulse 
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t = 20 / {q^/nQ) ^ where n is the number of phonons in the 
ground state of the given two-level system. Taking the 
variable to be of order 1, we need 0/2 ^ v/mnrj^ for 
each "vertical" two level system, and $7/2 ^ vrj^/n/m for 
each "diagonal" two level system, where m is the number 
of pulses applied to that given system while it is occupied. 
To make single inequalities for each two-level system, we 
take the maxima of mn = [{N + l)/2]^ and m/ y/n = N , 
where N is the maximum value of n in the state to be 
synthesized. 

We now define a "quality factor" q, and use this to 
calculate the laser Rabi frequencies VI / 2 = qAv/[{N + 
1)77]^ and ri/2 = qvq/N for the cases of the "vertical" 
and "diagonal" pulses respectively. These expressions 
replace the original inequalities described above, with q 
determining the "quality" of the laser pulse {i.e. the 
smaller the value of g, the greater the validity of the 
approximate Hamiltonians and (^), as appropriate). 
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FIG. 1. For all four plots: solid line, phase state with N = 8 
and = 2 radians, synthesized when q — 0.01; dotted line, 
same state synthesized when q = 0.1; dashed-dotted, phase 
state with A'^ = 32 and (jj — 2, synthesized when q — 0.01; 
dashed line, same state synthesized when q = 0.1. (a) Fi- 
delity of synthesis of the above states using a travelling wave 
laser configuration for a range of r;. (b) Fidelity of synthesis 
using a standing wave configuration, (c) Magnitude of the 
dimensionless quantity vt/2n required for the synthesis of the 
states over a range of t) (identical for both lasere configura- 
tions), (d) Same as (c), except using a logarithmic scale. 



close the synthesized state is to what it is supposed to be. 
Note that our calculations as to the circumstances when 
the use of the approximate Hamiltonians (||,||) is justified 
is restricted to the Lamb-Dicke limit, including our use 
of the quality factor q. We nevertheless continue to use 
q to define the laser Rabi frequencies up to 77 = 0.95 in 
Fig. 1^, for the sake of comparison. 

As indicated in Fig. |l], a completely analogous deriva- 
tion can be carried out for a standing wave laser configu- 
ration jsj . A "vertical" pulse is achieved when A = and 
the ion is at an antinode. A "diagonal" pulse is achieved 
when A = —v and the ion is at a node. The approxi- 
mate Hamiltonians are completely identical to the trav- 
elling wave case (^,^, except that e^"'(°+"^) is replaced 
by cos(77[a -I- a''']) and sin (77 [a -f a^]) for the "vertical" and 
"diagonal" pulses respectively. Thus, the constraints are 
somewhat different. In the case of the "vertical" pulses, 
when the cosine function is expanded to a power series 
there are obviously only even powers of 77(0 -I- d^). This 
means that when the "vertical" pulse is the pulse desired, 
the secondary pulse is not the "diagonal" pulse (which 
does not exist) but the pulse coupling \g, n) to |e, n — 2), 
which is much smaller. Similarly in the case of the "di- 
agonal" pulse, (sine expansion, and thus no odd powers 
of 7/(d -f a))) the secondary pulse is that coupling |(7,n) 
to |e,77 -|- 1), which is also small. Thus the constraints 
on the intensity of the laser are much less restrictive. 
We nevertheless use the travelling wave constraints for 
constructing the states shown in Fig. |l|, to show that sig- 
nificantly better results can be achieved using a standing 
wave configuration, where the circumstances are other- 
wise identical. Note also that we are generally used to 
thinking of the regime where the Lamb-Dicke limit can 
be assumed as being universally optimal. If we look at 
the plots in Fig. ^ however, we can see that higher values 
of F and lower values of ut are achieved for comparitively 
large 77. 
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The plots in Fig. |^ show the results of simu- 
lations creating Pegg-Barnett phase states — 

Y.n=oe"'^''"\^)I^NTl (where d^k = 2Trk/iN + l)) §. 
Below we will measure the Pegg-Barnett phase operator: 



N 



(10) 



fe=0 



and will therefore obviously need to know how to 
synthesize its eigenstates (the phase states described 
above). For these simulations we have used the ex- 
act Hamiltonian (^) for various values of q and 77. 
We define the fidelity F of the synthesis by F = 
|(desired state | synthesized state) p, as a measure of how 



FIG. 2. The filtering measurement process, (a) The ampli- 
tudes of all the ground states \g,n + 1) except for \g,0) are 
adiabatically transferred to the |e,n) states, (b) The laser is 
tuned to the decaying \g) <-> \r) transition. The presence of 
fluorescence will indicate that the ion was in the state \g,0) |^ . 
(c) If no photon is spontaneously emitted, then the original 
amplitude of the state \g,0) is projected out, and the popu- 
lation of the excited states \e,n) is adiabatically transferred 
back to the states \g,n + 1). 

The other tool we require in our measuring procedure 
is the ability to perform filtering measurements. First we 
shift all of the population in states \g, 77 4-1) to the states 
1 6,77) by adiabatic passage (where 77 > 0) so that 
all amplitudes with the exception of that for \g, 0) are 
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transferred unchanged (see Fig. ||). Next we use quan- 
tum jump techniques |5[| to determine whether jg, 0) is 
populated . Finally, if we do not detect population 
in state |5,0), we use adiabatic passage to restore the 
amplitudes in states |e,n) back to the states \g,n -I- 1). 
Adiabatic passage is required so that the population can 
be transferred independent of the Rabi frequencies. 
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FIG. 3. Plot showing the spread of the discrete posi- 
tion eigenstates l^fc) of Tig. x is in wavelengths/77, and 
y = |(a::|a;fc)|^. |a;o) is the leftmost peak, and Ixg) the right- 
most. 

We have simulated our measurement process, where 
one attempts to determine the phase and position ob- 
servables of particular states. To this end we used as 
bases: the eigenstates of the Pegg-Barnett phase opera- 
tor ( p^ ) 1^; and the eigenstates of the position operator 
xjv = ojv + ajy, for some truncated Hilbert space Hn+i- 
The position eigenvalues are = y for which the Hermite 
polynomial Hpf^i{y) = (see Fig. ||). For large N 
these zeros are separated by 2n/\/4N (asymptotic limit). 
The full Hamiltonian (^) was used in these simulations, 
and the results are shown in Fig. ^. It is especially in- 
teresting to consider the case of a superposition of two 
coherent states {\a) + \ — a)). As coherent states are 
quantum-mechanical representations of essentially clas- 
sical states, this can be considered a quantum superposi- 
tion of classical states, or "Schrodinger cat" state. Well 
separated peaks (in the position basis) can be obtained 
for values of a as low as 1.5. The overall "size" of such a 
state is very small however; typically the peaks are sepa- 
rated by less than a wavelength, so that they cannot be 
observed by just looking. Using our proposed measure- 
ment procedure, one would obtain in a single measure- 
ment one value of k that allows us to elucidate one peak 
or the other, as shown in Fig. ^ 

In summary, we have shown how to perform a single 
measurement for an arbitrary observable for the motion 
of a trapped ion. To our knowledge, this is the first ex- 
ample of an experimentally realistic system where arbi- 
trary measurements can be practically implemented. The 
method is based on being able to synthesize arbitrary 
states of motion in an ion trap, and performing filtering 
measurements on the state of the ion by quantum jump 
techniques. 
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FIG. 4. Measurement of the the phase (a,b) and posi- 
tion (c,d) observables (a,c) The ion is initially in the 
phase state with N — 8 and cj> = 2 radians; for these simu- 
lations q — 0.1 and rj = 0.5. (b,d) The ion is initially in the 
Schrodinger Cat state with A'^ = 32 and a = 1.5; for these 
simulations q = 0.1 and r] = 0.2. For each plot, black outline 
bars represent the results of simulation of the full Hamilto- 
nian, as given by Eq. ^ and grey solid bars represent ideal 
probabilities. 
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